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Abstract. Friction force acting on silver and nickel nanoparticles sheared on a 
graphene sheet in different lateral directions is investigated using classical molecular 
dynamics simulations. The results reveal the existence of frictional anisotropy for both 
metals. In most cases, the maximum value of the friction force is about two times 
larger than the minimum one. The form of dependencies of instantaneous values of 
the friction force components on the corresponding lateral components of the position 
of the centre of mass of the nanoislands strongly depends on sliding direction, varying 
between the sawtooth and the irregular one. A qualitative explanation of the results 
based on the "patch" model is proposed. 
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1. Introduction 

The ongoing endeavour to miniaturize various devices with moving components often 
meets an obstacle which hampers their reliable operability and which is related to friction 
at the atomic level [HEIGJIHEIEJIHIB]- A deep understanding of fundamental 
nanotribological processes is required in order to control friction and wear at the 
nanoscale and to successfully proceed with miniaturization tasks. In this context, 
manipulation of metal nanoparticles (NPs) adsorbed on an atomically-flat surface using 
the atomic force microscope (AFM) tip is worth mentioning. This method has been 
developed recently for quantitative investigation of friction of antimony NPs sliding on 
highly oriented pyrolytic graphite (HOPG) substrate [91 [TUl [HI [121 Q2]- Compared to 
traditional techniques it offers several advantages. First, the method provides much 
wider range of materials available for consideration. Second, with AFM it is possible to 
probe tribological properties of nanoobjects with sizes which are not accessible to friction 
force microscopy (FFM) and surface force apparatus. The manipulation experiments 
reveal the dual behaviour of the friction force Ff acting on Sb NPs under ultra-high 
vacuum conditions. Some NPs exhibit vanishing friction while others experience finite 
Ff which depends linearly on the contact area A [H]. Note that in the experiments 
typical diameter of a NP varies from about 50 nm up to a few 100 nm, and the contact 
area is deduced from lateral dimensions of the NP so the flat contact interface is assumed. 

From the theoretical side, several works pertaining to friction of metal NPs on a 
solid substrate can be found in literature [3 [TJ)J H5J HZJ Q3] . Seminumerical model [TJ] is 
aimed at calculating the power necessary to depin a NP on a surface with an AFM tip 
in tapping mode. However, a large amount of assumptions, such as, two-dimensionality 
of the system, the rigidity of the NP and the substrate, quasistatic motion and the 
absence of the Ff measurements make it inappropriate for comparisons to the kinetic 
friction experiments. Molecular dynamics (MD) simulations of the ballistic friction [7] 
consider sliding of a golden NP on the HOPG substrate with high velocities and do not 
measure Ff. Therefore, this model cannot describe the manipulation experiments. Ab 
initio calculations of friction of a Sb4 molecule on HOPG surface [33] make an attempt 
to explain the dual behaviour of Ff in the experiments. However, it is questionable 
whether the results for one Sb molecule can be transferred to the relatively large NPs 
comprising hundreds of thousands and more atoms. In addition, sliding in a static 
manner [TB] makes it questionable whether it is valid to compare the calculations [TB] 
with the experiments, where NPs move with velocities of the order of /im/s. 

MD simulations of friction of silver and nickel [T7] as well as golden and copper [T5] 
nanoislands adsorbed on a graphene sheet provide a model which is relatively close to 
the manipulation experiments. Although the typical lateral dimension of the largest 
NP is about 10 nm which is 5 times smaller than in the quantitative tribological 
experiments, such NPs are used in other manipulation trials [T9J [201 EI]- Further 
enhancement of sensitivity of the equipment will give the opportunity to perform 
quantitative measurements for such nanoislands and to compare with computational 
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results. Another deviation of the model p2] from the experiments is the use of different 
materials, viz graphene as the substrate and Ag or Ni for the NPs. Nevertheless, such an 
approach should reproduce the general behaviour of the experimental systems, because 
their features which are essential to tribological properties should not differ significantly 
from our model. For example, the structural contribution of graphene to friction of 
nanoislands should be close to that of HOPG because these materials have identical 
lattices. Friction and other tribological properties of graphene itself have attracted 
significant attention in the last years (221 E3J EU |25j |26] and our simulations may 
provide valuable information about it. Additionally, the development of new methods of 
production of Ag [27J ES] , Ni [25] and Pt [30J NPs with sizes comparable to our system 
provide the possibility for experiments which directly correspond to our model. 

The main aim of the simulations in Refs [T7J US] was to define the dependence 
of the friction force acting on NPs on contact area. However, many other unresolved 
questions have remained. One of them is the dependence of the friction of nanoislands 
on the sliding direction relatively to the honeycomb lattice of graphene. The 
existence of the frictional anisotropy in different systems has attracted a significant 
attention recently [3H E21 E3J EH E5] because it allows one to separate the dominating 
contributions to friction from a number of other factors. It can also help clarify the 
origin of the superlubricity phenomenon [331 EE]- According to the experiments [31 J and 
simulations [321 E3J El], frictional anisotropy usually exists in small systems where the 
contribution of the atomic structure of the interface cannot be neglected. However, it 
is not necessary for both contacting surfaces to have periodic crystal lattice in order 
to observe the phenomenon as was shown for quasicrystals [3H EE] and lubricated 
systems [31]. One should note that the results mentioned above were obtained for 
kinetic friction, while static friction in the one-atom Prandtl-Tomlinson model for the 
hexagonal lattice of the substrate is isotropic [33J . 

According to the manipulation experiments, Sb nanoislands do not manifest 
frictional anisotropy on HOPG surface [H]. Additionally, no locked state of NPs was 
found during their rotation. These results may be caused by the following two factors. 
First, relatively large sizes of nanoislands with hundreds of thousands of atoms in contact 
may lead to washing out of the anisotropy effect. Second, the nearest-neighbour distance 
in Sb is much larger than in graphene and may prevent the interface from forming 
locked regions. In order to clarify these questions, we have carried out MD simulations 
described in the present paper, where the direction dependence of the kinetic friction of 
silver or nickel nanoisland sliding on graphene is investigated. 

2. Model 

The details of the model can be found elsewhere [UJ E7J. Here we point out only its 
main peculiarities. A graphene sample has zigzag and armchair edges parallel to the x- 
and ^-directions, respectively (see figure [T]). To maintain the constant position of the 
graphene in the ^-direction, boundary carbon atoms along the perimeter of the layer 
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are kept rigid in the simulations. Ni and Ag NPs containing 10000, 20000, 25000 and 
29000 atoms are considered. For each size of the NP the unique x x y dimensions of the 
substrate are used, varying from 20.7x17.9 nm to 37.4x32.4 nm. The total number of 
atoms in the simulations is in the range from 24112 to 75208. 




Figure 1. Top view of the system with a Ni nanoisland containing 29000 atoms (left). 
Schematic sketch of the forces acting on the NP during its shear (right). (All snapshots 
in this work were produced with VMD software 38J). 

Graphene atoms interact with each other via the harmonic potential [3S]. Forces 
between metal atoms are derived from the alloy form of the embedded atom method 
potential [30]. The pairwise 6-12 Lennard- Jones potential with the parameters from [T7] 
describes the metal-carbon interactions. The equations of motion are integrated using 
the leapfrog method [311 132] with a time step At = 0.2 fs. Calculations were performed 
on GPUs NVIDIA GeForce GTX 460 and GeForce GTX 480. 

A general flow of a typical simulation can be deduced from the figure [2] A 
nanoisland is obtained by the procedure imitating the dewetting of thin metallic films by 
thermal treatments [T71 l29"j 14*3] . In this approach, a thin slab of metal atoms packed into 
an ideal fee lattice is placed initially on the graphene. After the start of the calculations, 
it self-organizes into more compact nanoisland through the exothermal process. When 
the lateral sizes Lx and Ly of the NP diminish enough and the NP gains a desired 
symmetric elliptical shape, we apply the Berendsen thermostat [42J both to metal and 
the graphene during a suitable time interval to cool the system down to the temperature 
of about 300 K. After that the thermostat is decoupled from the metal atoms and is 
applied only to the graphene to dissipate the heat generated during the shear of the NP. 
The sizes Lx and Ly are computed as the difference between the coordinates of metal 
atoms with maximum and minimum values along the corresponding direction. Values of 
L x x L Y are about 6.7x6.3, 7.6x8.0, 8.4x8.5, 8.9x8.9 nm for Ni NPs containing 10000, 
20000, 25000, 29000 atoms, respectively. For corresponding Ag NPs these quantities are 
7.1x7.2, 8.6x9.1, 9.0x10.0, 10.0x10.2 nm, respectively. 

After the nanoisland has formed, it is sheared by the force F a imitating pushing 
with the AFM tip in the manipulation experiments. Shear force is applied to every 
atom located in the hatched sector in the right figure [T] along the specified direction 
defined by the angle 9 (counted from the x direction). The sector is formed by the lines 
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Figure 2. Time dependencies of quantities measured in the simulation of Ni NP 
containing 10000 atoms sheared in the direction corresponding to 9 = 30°. Here, Fx, 
Fy - lateral components of the substrate force; Lx _ lateral size of the NP in the 
x-direction; Xcm , Vx - x-componcnt of the position and the velocity of the centre of 
mass of the NP, respectively; T - temperature of the system. In the topmost plot the 
horizontal line is the applied force F a and the second (irregular) one is Fx- 



1 and 2, and the angles a — /3 — 30° provide approximately the same number of atoms 
in each half of the sector (symmetrically to the shear direction). Such a geometry gives 
zero torque relative to the NPs centre of mass (CM) in the xy-pl&ne, thus preventing 
the NP from rotating in this plane. A metal atom belongs to the hatched area if its 
coordinates satisfy the conditions following from the equations of the lines 1 and 2: 



y < xtan (6 — a) 
y > x tan (9 + 0) 



(1) 



At the beginning of a simulation run, x- and ^-components of the applied force 
are incremented in steps of 0.0001 pN until the x-component of the velocity of the CM 
Vx of the NP reaches the value of 3.55 m/s. Then the applied force acting on each 
atom in the sector remains constant, and the simulations are held with constant total 
applied force F a . In the computer experiments both lateral components Fx and Fy of 
the substrate force and the coordinates of the CM of the NP Xcm, ^cm are measured. 
Each force component is defined as the sum of the corresponding components of forces 
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acting on metal atoms from the graphene. The total friction force Ff = J (Fx) 2 + (Fy) 2 
is defined by the time-averaged values of the lateral components of the substrate force. 
We consider the following values of 9: 0°, 22.5°, 30°, 45°, 60°, 67.5° and (for some NPs) 
89°. One should note that MD simulations do not allow one to investigate the low 
NP sliding velocities (of order ~ /im/s) used in the AFM experiments. However, the 
simulation results obtained in Refs. [TSJ HI] suggest a very weak dependence on velocity 
in the studied range. 

3. Results 

Let us analyze the dependencies of the total friction force Ff on the angle 9 displayed 
in figures [3j [4], [5j Markers in the figures correspond to the calculated values and the 
solid lines are cubic splines shown to follow the eye. Generally, more or less pronounced 
changes of Ff with 9 are observed for all considered NPs. For most cases, the maximum 
value of the friction force is about two times larger than the minimum one. There is 
a correlation in variations of Ff with 9 for Ag NPs containing 10000 and 20000 atoms 
(figure [3j left part), with local maxima in the angle range between about 18° and 30°. 
However, for the corresponding Ni NPs the picture is different with considerable data 
scatter (cf. figure |3j right part). Such a result can be attributed to the fact that the 
nearest-neighbor distance c?Ni = 0.249 nm in Ni is very close to the graphene lattice 
constant a = 0.246 nm. This may lead to a high sensitivity of Ff to small variations 
of number of metal atoms in contact as there is a high probability of the formation of 
locally commensurable interface leading to high friction [T71 [18]. The procedure used 
in our model for production of nanoislands did not allow us to precisely control the 
shape and the area of contact of NPs, and the random factor is always present. So NPs 
containing the same number of atoms do not have the identical contact morphologies in 
different simulation runs. This fact for small Ni NPs may lead to considerable deviations 
in the friction force. 

Data shown in figure [4] also indicates the existence of variations of F{ with the 
direction of sliding for larger NPs. The maximum and the minimum values of F{ differ 
by about factor of 2. This value is typical for many systems, but it is smaller than the 
one found for quasicrystal substrate, where it is equal to 8 [31J. Data scatter did not 
allow us to certainly extract the empirical dependence Ff(9). Considerable variations 
of Ff with misfit angles of substrates were also observed in the lubricated systems j3l] 
especially for the thinnest, one atom thick lubricant film. In this case, the origin of 
variations is the commensurability of the lubricant-substrate interface due to solidlike 
or liquidlike state of the film. In our model the reason of variations of friction force is 
different, as the contact interface of the NP does not have long-range ordering of atoms 
(see below). Data scatter may be ascribed to slightly different contact morphologies for 
NPs containing the same number of atoms in different simulation runs, which may lead 
to random occurrences of uncontrollable locally commensurable contacts. For Ag and 
Ni nanoislands composed of 29000 atoms qualitatively similar trend of F{(9) dependence 



Frictional anisotropy of metal nanoparticles adsorbed on graphene 



7 



0.28 



0.24 - 



0.16 




0.12 - 



0.08 



0.04 




100 



Figure 3. The total friction force acting on Ag and Ni nanoislands containing 10000 
and 20000 atoms for different sliding directions. 



can be observed (figure [5]) with the minimum between 30° and 45°. 





Figure 4. Direction dependencies of the total friction force acting on nanoislands 
containing 25000 and 29000 atoms. 



Let us analyze the dependencies of the lateral (substrate) force components Fx and 
F Y on the corresponding components X C m, Ycm of the CM of NPs shown in figure [6] 
and figure [7| In the previous study [T7j it was found that Fx(Xcm) has an irregular 
shape for Ag nanoislands sheared along the direction corresponding to 6 = 0, while 
this dependence is sawtooth or has sawtooth regions for most Ni NPs. In the present 
work we show that both Fx and Fy components can be sawtooth or have such regions 
for both metals when 6^0°. This can already be seen from the time dependencies in 
figure [2j where both Fx and Fy are sawtooth for Ni nanoisland containing 10000 atoms 
and 9 = 30°. Most dependencies are sawtooth or have such regions for Ni NPs, while 
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Figure 5. Total friction force vs angle 9 for Ag and Ni NPs containing 29000 atoms. 

for Ag such a dependence is an exception rather than a rule. 

As can be seen from figure [6]a, b, shearing Ag NPs with 29000 or 20000 atoms in 
the direction corresponding to 9 = 22.5° leads to sawtooth components of the friction 
force, while other components (not shown here) are not clearly sawtooth. The distance 
between spikes fluctuates around the value of about 0.5 nm. Previously [T7] it was 
shown that the period of spikes for Ni was close to the value of the lattice constant a 
of graphene. Deviation of the shear direction of Ag nanoisland with 29000 atoms from 
9 = 22.5° to 9 = 30° leads to the qualitative change of the dependencies (figure [6]c, d) 
from sawtooth to the irregular one. 

In the case of Ni nanoislands the dependence of at least one of the components of 
the friction force is sawtooth or has sawtooth regions for the majority of calculations. In 
some simulation runs both components are sawtooth as is the case for Ni NPs with 10000 
and 25000 and 9 = 30°, while only one component is not irregular in others. In general, 
the period of spikes differs from the Ag NPs having the value of about 0.25 nm close to 
a[T7j. 

Our results clearly indicate that the friction force changes with the sliding direction. 
To explain this fact, several factors can be considered. The first one is the structure 
of the contact interface. We calculated radial distribution functions (RDFs) to look 
at the structure of nanoislands. For all NPs RDFs are analogous to the obtained 
previously p21 [18] (they are not shown here) and have the first spike indicating the 
short-range order in NPs. The location of the spike is close to the nearest-neighbor 
distance 0.249 nm and 0.288 nm for Ni and Ag, respectively. Farther spikes are 
smeared, indicating the absence of the long-range order and suggesting that the bulk 
of NPs is either amorphous or polycrystalline. Visual inspection of the contact surface 
of nanoislands confirms this conclusion and suggests that the NP interface contains 
domains with crystal structure and some disordered regions. 
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Figure 6. Examples of sawtooth (a, b) and irregular (c, d) dependencies of the 
substrate force components for Ag NPs. 



The second factor is the influence of the energy dissipation mechanisms. According 
to the experiments [31] , long-range order of both surfaces is not required to observe the 
frictional anisotropy, and it can occur when one of the sliding surfaces is amorphous. 
This was shown for TiN AFM tip passivated with a molecular layer of hexadecane 
thiol and sliding over the Al-Ni-Co single quasicrystal surface. The authors found that 
the friction along the aperiodic direction of the surface is eight times smaller than the 
force along the periodic direction. They explain such a behaviour by energy dissipation 
through electronic and phononic mechanisms, which are known to be highly anisotropic 
in the considered quasicrystals. The authors rule out the structural contribution due 
to an anisotropic response of the hexadecane thiol molecules that coat the tip, because 
the difference between vertical corrugations along periodic and aperiodic directions is 
small compared to the size of the alkyl chains. They also state that incommensurability 
between the probe and the surface is unprobable in any scan direction because the 
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Figure 7. Sawtooth (a, b) and irregular (c, d) dependencies of the substrate force 
components for Ni NPs. 



TiN tip is amorphous and is covered by thiol molecules. In contrast to the reasoning 
presented above, the simulations suggest that the commensurability factor can play an 
important role in similar systems [32j Ef5] . 

In our case the mentioned factors can also be considered in order to explain 
the obtained results. However, unlike the Ref. [31], the phononic and electronic 
contributions cannot lead to anisotropic effects because in our model they are taken 
into account through the Berendsen thermostat. The latter acts on all carbon atoms, 
so the energy is dissipated isotropically. Thus, the structural contribution is the most 
significant one in our model. 

A "patch" model can be employed to qualitatively explain the results. It was 
introduced in Ref. [18J to find the reasons for the sawtooth dependencies of the friction 
force for NPs sliding along the x-direction. The key point of the model is the assumption 
that the friction force is mainly defined by the ordered regions of atoms on the surface 
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Figure 8. Illustration of contact "patches" for Ni NP containing 25000 atoms and 
sheared at angle 60° at time moments of 106 and 132 ps. Patches are outlined by 
contours. 

of the NP. These ordered domains can form commensurate regions with the substrate 
(see figure [8]) which give the major contribution to the friction. In order to get high 
friction and sawtooth substrate force dependencies, the spots should be easily formed 
and act coherently during the sliding of the NP. This was observed for Ni [T7j and 
Cu [18] where the formation of the "patches" is facilitated by a good fit of the metal 
nearest-neighbor distance to the graphene lattice. On the other hand, for Ag and Au 
the formation of patches is more difficult and hence the sawtooth dependencies were 
observed only as exception at low sliding velocities. The anisotropy of the friction force 
can also be attributed to the dependence of the number and the topography of patches 
on the sliding direction. Those directions where the number of patches is larger are 
characterized by higher friction. Additionally, in some directions the patches do not 
act coherently leading to irregular forces, while in other directions the coherent action 
leads to sawtooth curves. This reasoning is consistent with our simulations, where the 
force curves are either irregular or sawtooth for Ag depending on the sliding direction 
of the NP. Quantification of the "patch" model and its application to the corresponding 
systems will be presented in future studies. 

4. Conclusion 

In conclusion, classical MD simulations have revealed the dependence of the friction force 
on the shearing direction for Ni and Ag nanoislands. It was shown that the anisotropy 
effect can lead to the difference in friction of a factor of 2. The dependencies of the 
components of the substrate force on the corresponding lateral components of the CM 
coordinates of the NP are also changed with the shear direction. The possibility of the 
sawtooth dependency was shown for Ag NPs. In this case the frequency of spikes differs 
from the one of the sawtooth friction curves for Ni. A qualitative "patch" model was used 
as a first effort to explain the friction anisotropy phenomenon. Further advancement of 
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the model to the quantitative level and its application to different systems will be done 
in future studies. 
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